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Abstract 

We present the Mediatrix filamentation method, a novell iterative 
procedure that decomposes elongated objects in filaments along their 
main direction over their intensity peak. From this decomposition, 
(-H I the method measures the object's length and thickness. This tech- 

04 nique is applied in preliminary tests to arc-shaped objects (simulated 

gravitational arcs) to recover their curvature center. 
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1. INTRODUCTION 



Shape analysis and detection are fundamental issues in image processing. 
For curved shapes, it is of interest to define a curvature, with an associ- 
ated center and radius. We may find this kind of shape, among others, in 
astronomical images. The quintessential example in this context is that of 
gravitational arcs. These objects are the result of the Strong Lensing ef- 
fect P, |2] which occurs when a distant source object is aligned with some 
intervening distribution of matter in the trajectory of light towards us. This 
matter distribution distorts the space-time, acting like a lens. The image of 
astronomical objects is distorted and magnified, forming the curved shapes 
mentioned above. We have developed a novell technique to decompose images 
of elongated objects in a set of filaments and measure their length, width, 
and assign a curvature center. 

This paper is organized as follows: In section we introduce the Me- 
diatrix Filamentation method explaining its basic elements and procedure. 
Section SL presents measurements derived from the Mediatrix method. In 



^E-mail: debom@cbpf.br 



1 



section 4, we apply the method on simulated arcs and discuss the results. In 
section 5^, we present our concluding remarks and future perspectives. 



2. MEDIATRIX FILAMENTATION METHOD 

In the following we will describe the method to decompose elongated 
objects into a set of line segments. For concreteness we will consider the 
object to be composed by a set of pixels with given intensity, as for a digital 
image. However, in principle, the method can be applied to any intensity 
distribution, even if not pixelated. The only requirement is that the object 
should have a clearly defined boundary, in other words, the intensity must 
be zero outside the object. 

The Mediatrix method was originally designed to search for gravitational 
arcs. It was inspired on a basic geometrical property of the perpendicular 
bisector of pairs of points on a circle. Namelly that these lines, for any set 
of pairs of points, intersect at the circle center. Therefore if an elongated 
object can be decomposed into a set of points along its longer direction, and 
if this object has a shape close to an arc segment the perpendicular bisectors 
of pairs of these points will intersect in nearby points (i.e. close to the center 
of curvature). 

It turns out, however, that this method can be used to assign segments 
along the longer direction of any elongated object, i.e. to "filament" the 
object, or to determine its "spine", regardless of the presence of curvature. 
The key procedure to filament the object is to recursively obtain the perpen- 
dicular bisector of pairs of points on the image. 

Given the points Pi = {xi,yi) and P2 = (3:2, 2/2)5 the perpendicular bi- 
sector is a straight line y = mx + h perpendicular to P1P2 that crosses the 
middle point of this segment and whose coefficients given by: 



The Mediatrix method is a recursive method that operates on several 
iteration steps. Each step is a new Mediatrix level and the method can be 
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performed up to an arbitrary level n. In the following, we describe the first 
few levels as an example. 

In the first step we determine the extreme points, Ei and E2 of the object. 
Several methods have been considered to determine the extreme points of an 
object (see, e.g., ref [3]). Here we use the "farthest-of-farthest" method, 
by which Ei is defined as the most distant point from a reference point on 
the object (e.g., the brightest pixel on the image or its geometrical center), 
whereas E2 is defined as the pixel on the object farthest from Ei. Next, 
the first perpendicular bisector of these two points is calculated. The first 
Mediatrix point is defined as the brightest pixel of the object along the 
perpendicular bisector. In pratice, we take the brightest pixel located at a 
distance d < aAp from the perpendicular bisector, where Ap is the pixel 
scale and a is a free parameter choosen as a = The first Mediatrix Point 
is shown in Fig. l(A)[T]for an arc-shaped object (more specifically, an 
ArcEllipse [1]). 

In the second step the perpendicular bisector is now calculated with re- 
spect to El, and M^, E2. These two perpendicular bisectors now define 
two other Mediatrix Points: and M| using the same criteria we used to 
define (Fig IB). The upper index refers to the iteration level and the 
lower index is a label to identify the points. Proceeding to the third step, 
presented in Fig 1(C), we have the set of Mediatrix Points M^, Mf, M| and 
the two extremes Ei and E2. Those points are used to define new Media- 
trix Points Mf calculated by picking the highest intensity pixel near to the 
perpendicular bisector between two neighboring points. The algorithm con- 
tinues defining new Mediatrix Points M- , corresponding to the i-est point in 
the j-est iteration level, in higher iteration levels until reaching a specified 
final step n. In Fig 1(D), we present the last step for n = 3 (as in the previ- 
ous panels some points were omitted not to crowd the figure). The collection 
of Mediatrix Points together with the two extreme points are then named 
keydots. From the keydots, the object is decomposed in = 2" segments or 
filaments. Each segment connects a keydot to its neighbors. The algorithm 
outputs a set of vectors nj, where j varies from 1 to A^. Those vectors are 
perpendicular to the segment which connects a keydot to its neighbor with 
origin in the middle point of its segment and norm equal to the length of this 
segment. This is shown in Fig. 1(D) for n-r, where |n7| = |M4 Ml]. 
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3. MEASUREMENTS DERIVED FROM THE MEDIATRIX 
METHOD 



Using the outputs of the proposed algorithm, the object length, L, is 
defined as: 



N 



L = j:\n,\. (3) 



The points Fj and Fj, represented in Fig. 1(D) for j = 7, are defined as 
the two extreme pixels from the set of points along the perpendicular bisector 
associated to Uj. Using those points, it is possible to measure the average 
thickness W: 

w=^hFhF^ (4) 

For arcs constructed from circle segments, all perpendicular bisectors in- 
tercept at the center of curvature. Therefore, we may define as the object 
center of curvature the average of the points Cik = {xik, Uik) generated by the 
intersection of the lines in the nj and directions. The center coodinates 
Ca = {xa, Va) are given by: 
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i=l kf^i 



Va = 2^ 2^ Vik (6) 

■ i=l fc^i 

An alternative way to define a curvature center is using the median inter- 
cept point Cm = {xm,ym), wherc Xm is defined as the median of the values 
of Xik and i/m is given by the median of i/ik- 

From these two "center of curvature" definitions we may assign a cur- 
vature radius Ra and Rm, associated to Cm and Ca, as the distance from 
the first Mediatrix filamentation point to Cm or Ca respectively, i.e., 
Ra = IWC^I and Rm = \M^\. 
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The Mediatrix filament at ion method was implemented in the Python pro- 
gramming language [5]. The code was developed as part of the SLtools li- 
brarj|§. 

4. APPLICATION TO SIMULATED GRAVITATIONAL ARCS 



As an example of application of the Mediatrix filamentation method, 
we consider 3 arc-shaped objects produced through the gravitational lensing 
effect (see Fig. 2|2]). These arcs were generated using the AddArcs [6] pipeline 
to simulate gravitational arcs. 

AddArcs is a code that simulates realistic gravitational arcs using the 
galaxy cluster abundance (i.e. the lenses) provided by cosmological simula- 
tions and background galaxies (the sources) with morphological parameters 
and redshift distribution obtained from the Hubble Ultra Deep Field Survey. 
The source brightness distribution is modeled by Sersic profiles [7], which in 
principle extend to infinity, with elliptical isophotes. Given the input models 
for the source and the lens, AddArcs controls the software gravlens [5] to 
perform the gravitational lensing calculations. The segmentation of the sim- 
ulated image (i.e. the definition of borders) is carried out with the SExtractor 
[9] software, in a procedure similar to the one described in [1] 

To perform the Mediatrix filamentation we used the function 
run_mediatrix_decomposition_on_stamp from the mediatrix_decomposition 
module in SLtools, which was developed to apply the method to an image 
matrix with one isolated arc. The points {xik,yik) are shown for the tested 
arcs in Fig. 3. In order to avoid lines that do not intercept or are almost par- 
allel to each other, the algorithm ignores combinations of Ui and which 
are almost in the same direction, more specifically, we discard pairs with 
\tan6i — tan6k\ < 10~^, where 6i and 6k are the angles determined by Ui and 
fik directions with respect to the x axis. The center of the mass distribution 
that generated the the lensing effect, the lens center Ci = {xi,yi) is given 
by AddArcs. Although from gravitational lensing theory there is no need 
for the curvature center and the center of mass to coincide, they are usually 
close. We define the arc radius R as the distance from the first Mediatrix 

^This pipeline is part of SLtools, a library for image processing, catalog manipulation, 
and strong lensing applications [lOj . and is available at http://che.cbpf.br/sltools/. 
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u 


R 


L (pixels) 


L/W 


Arc A 


4 0.036 


0.062 


294.7 


5.6 


8 0.153 


0.126 


295.3 


5.6 


16 0.144 


0.092 


295.7 


5.6 


32 1.09 


0.230 


296.8 


5.7 


Arc C 


4 0.100 


0.096 


362.4 


12.6 


8 0.211 


0.68 


363.1 


12.7 


16 0.166 


0.170 


363.3 


12.7 


32 0.307 


0.168 


363.9 


12.7 


Arc C 


4 0.101 


0.097 


192.9 


5.5 


8 0.051 


0.093 


193.1 


5.5 


16 0.085 


0.126 


193.5 


5.6 


32 0.198 


0.209 


194.7 


5.6 



Table 1: Results of Mediatrix measurements for 3 simulated arcs A, B and 
C. The first column is the number of segments the arc was divided into. The 
second and third are the ratio between the distance from the lens center 
to the center of curvature — calculated using the average ARa and median 
ARfn respectively — and the curvature radius R. The fourth column is the 
arc length and the fifth is the L/W ratio. 



filamentation point to C;, so i? = |M^C/|. The distance between Ci and 
the measured centers of curvature Cm or Ca is AR^ and ARa respectively. 
We present the results for Ra, Rm, discrepancy in R using R^, L, L/W 
and curvature R/L'm Tabled] for the 3 arcs. 

4.CONCLUSIONS AND FUTURE APPLICATIONS 

The Mediatrix filamentation method is a technique that decomposes and 
measures elongated shapes, and may be used in the detection of shapes such 
as gravitational arcs. The technique can also define an expected center for 
curved objects. Previous definitions of the curvature center in the gravita- 
tional arc characterization problem assumed that the arcs are circular [TT] : 
with the Mediatrix method there is no need to use this hypothesis. 
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The rij vectors can be used as input for other morphological estimators for 
arc-shaped objects. For example, in ref. [1] the Mediatrix output is used as a 
starting point for a method that fits the brightness distribution of arcs using 
analytical templates to derive structural parameters of the arcs. Another 
appication is to use parameters derived from the Mediatrix decomposition as 
input for a neural network to classify objects fn\ . 

In the current case, the extreme points were used as the first Mediatrix 
step but other definitions for the first step can be used, e.g. object edges, in 
order to filament more complex objects not necessarily with a single preferred 
direction. 
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(A) (B) 

E2 E2 




Step n: Define new Mediatrix Final Step:The algorithm assigns a 

Points between each pair of set of vectors perpendicular to the 

neighbors Mediatrix Points from segments that connect two 

all previous Steps. neighboring keydots. 



Figure 1: Steps in the Mediatrix Filamentation method. After n iterations, 
the method determines a set of 2" points defined by the maximum of in- 
tensity along the 2" perpendicular bisector and the 2" vectors perpendicular 
to neighboring points with magnitude given by the distance between these 
points. For clarity, only some points afe shown on the figure, which illustrates 
the steps for n = 3. 



Figure 2: Arcs A, B, and C respectively used as input for the Mediatrix 
method. 
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